STATUS OF CENTER DOMINANCE IN VARIOUS CENTER GAUGES 
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Abstract. We review arguments for center dominance in center gauges where 
vortex locations are correctly identified. We introduce an appealing interpretation 
of the maximal center gauge, discuss problems with Gribov copies, and a cure 
to the problems through the direct Laplacian center gauge. We study correlations 
between direct and indirect Laplacian center gauges. 



1. Why Center Dominance? 

The aim of most lattice studies of the confinement mechanism is to extract from 
lattice link variables the most relevant parts for the infrared dynamics. The con- 
cept of (some kind of) dominance seems a necessary, though not sufficient, condi- 
tion for success. If we extract the (would-be) relevant parts of links and compute 
physical quantities related to confinement (e.g. the string tension) we expect to 
reproduce their behavior in the full theory. Were it not the case, one could hardly 
claim to have achieved the goal. 
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It was observed both in SU(2) [1,2] and (less convincingly) in SU(3) lattice 
gauge theory [3] that the string tension obtained from center-projected configu- 
rations in maximal center gauge (MCG) agrees remarkably well with the asymp- 
totic string tension of the full theory. This phenomenon of center dominance has 
led to the recent revival of interest in the center-vortex picture of color confine- 
ment [4-6]. One can easily formulate an argument why center dominance should 
occur if center vortices are correctly identified [2] : 

Vortices are created by discontinuous gauge transformations. Let a closed 
loop C, parametrized by x^{x), 1 € [0, 1], encircle n vortices. At the point of dis- 
continuity (in SU(2)): 

g{xm = {-irg{x{\)). (1) 

The corresponding vector potential in the neighborhood of C can be decomposed 
as 

aJ") (x) = g-i5A W {x)g + ig-%g. (2) 

The g^^d/jg term is dropped at the discontinuity. Then, the value of the Wilson 
loop is 

W„{C) = (Trexp[/^ Ja^aJ,"']) = (-l)"(Trexp[/^ J^SaJ,"^]). (3) 

In the region of the loop C, the vortex background looks locally Uke a gauge 
transformation. If all other fluctuations SA^"'' are basically short-range, then they 
should be obhvious, in the neighborhood of the loop C, to the presence or absence 
of vortices in the middle of the loop. In that case: 

(Trexp[j^^/y5Ai"^) ^ {TTexp[i ^ dx^'dA^]) (4) 

for sufficiently large loops, and therefore 

Wn{C)/Wo{C) — ^(-1)" or W(C)!^Z(C)x(Trexp[j^jy6A{,°^]). (5) 

Here W{C) is the expectation value of the full Wilson loop, and Z{C) the expec- 
tation value of the loop constructed from center elements alone. 

It is clear from the above argument that one gets center dominance, i.e. the 
same string tension from W(C) and Z(C), under four intertwined assumptions: 

1 . Vortices are the confinement mechanism. 

2. Vortices are correctly identified. 

3. Short-range fluctuations with/without vortices look similar. 

4. No area law arises from the last factor in Eq. (5). 

Is there a necessity to fix any gauge? The original vortex idea was formulated 
without a reference to a particular gauge, and in fact a kind of center dominance 
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exists even without gauge fixing, as was shown in [7]. However, this holds for 
any distances, not only for large ones, vortices defined without gauge fixing do 
not fulfill simple expectations and do not scale according to the renormalization 
group, and thus the phenomenon hardly bears any information on the confinement 
mechanism. Gauge fixing appears of special importance for correct identification 
of vortices. 



2. How to Identify Center Vortices? 

The procedure, proposed in [1, 2], consists of three steps: 

1. Fix thermahzed SU(2) lattice configurations to direct maximal center (or ad- 
joint Landau) gauge by maximizing the expression: 

£|Tr[?7^(x)]|' orequivalently Y.^r[U^{x)] . (6) 

2. Make center projection by replacing: 

U^{x) ^ Z^(x) = sign Tr[U^ix)] . (7) 

3. Finally, identify excitations {P -vortices) of the resulting Z2 lattice configura- 
tions. 

A whole series of results, obtained by our and other groups, indicates that 
center vortices defined in MCG play a crucial role in the confinement mechanism. 
This includes, besides center dominance, the following: 

1. P- vortices locate center vortices in full lattice configurations [2]. 

2. P-vortices locate physical objects, their density scales according to the renor- 
malization group [8]. 

3. Creutz ratios computed from center-projected Wilson loops are almost con- 
stant starting from shortest distances; the Coulomb contribution was effec- 
tively eUminated (precocious linearity) [2]. 

4. Center vortices are correlated not only with confinement, but with chiral sym- 
metry breaking and non-trivial topology as well [9]. 

5. Deconfinement can be understood as a center vortex percolation transition 
[10,11]. 

Other gauges work as well; general conditions a suitable gauge has to fulfill 
were formulated in [12]. Here we would just like to briefly summarize an interest- 
ing insight into the meaning of MCG fixing and center projection, due to [13, 14]. 

3. Best-fit Interpretation of MCG 

Running a MC simulation, one can ask for the pure gauge configuration closest, 
in configuration space, to a given lattice gauge field: 

U^{x) - g{x)g\x+n)=ut-'\x). (8) 
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It is easy to show that finding the optimal g{x) is equivalent to the problem of 
fixing to the Landau gauge. 

Let us now allow for Z2 dislocations in the gauge transformation, i.e. fit the 
lattice configuration by a thin center vortex configuration: 

f/7'-(x) = g{x)Z,{x)gHx+fi), Z^{x) = ±1 (9) 

UJ,°''{x) becomes a continuous pure gauge in the adjoint representation, blind to 
the Zu(x) factor. One can make the fit in two steps: 

1. Determine g{x) up to aZa transformation by minimizing the square distance 
d\ between JJ^i and UJ,'"' in the adjoint representation, which is easily seen to be 
equivalent to fixing to direct MCG. 

2. Find Z^{x) by minimizing: 

Tr{ [U,{x)-g{x)Z,{x)g\x+ti)\ [u;{x)-g{x + m,{x)8Hx)\], (10) 

which requires the center projection prescription, Eq. (7). 

Summarizing, the procedure of direct MCG fixing + center projection repre- 
sents the best fit of a lattice configuration by a set of thin center vortices. 

4. Why Does MCG Sometimes FaO to Find Vortices? 

MCG fixing suffers from the Gribov copy problem. The iterative gauge-fixing 
procedure converges to a local maximum which will be slightly different for every 
gauge copy of a given lattice configuration. 

At the first sight, the problem seemed quite innocuous: We observed in [2] 
that vortex locations in random copies of a given configuration were strongly 
correlated. However, the successes of the approach were seriously questioned. 
Bomyakov et al. [15] showed that using the method of simulated anneaUng in- 
stead of our usual (over-)relaxation, one could find better (local) MCG maxima, 
but the center-projected string tension was only about 2/3 of the full one. 

The best-fit interpretation of the previous section provides us with a clue to 
the origin of this problem. It is clear that U'^°''{x) is a bad fit to Ufj{x) at links 
belonging to thin vortices (i.e. to the P-plaquettes formed from Z/jix)). We recall 
that a plaquette p is a P-plaquette iff Z{p) = — I (where Z{C) denotes the product 
of Zij{x) around the contour C) and that P-plaquettes belong to P- vortices. Let us 
write the gauge transformed configuration as 

W^{x) = Z^{x) e^" W , Tr e^'-W > 0. (11) 
At large P values ^Tr[Up] = 1 - 0(1/|3), and equals to 

(Zp) iTrn^^^W = p_p,,^,,,,, (-1) X iTrn^^^W. (12) 
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The last equation implies that at least at one link belonging to the P-plaquette 
A^(x) cannot be small, therefore ^Ujj{x) must strongly deviate from the center 
element. 

The above argument shows that the quest for the global maximum may not 
always be the best strategy; one should rather try to exclude contributions from 
P-plaquettes where the fit is inevitably bad [14], or modify the gauge fixing pro- 
cedure to soften the fit at vortex cores. 

5. A Cure for the Disease: Direct Laplacian Center Gauge 

We have recently proposed to overcome the Gribov problem using the direct 
Laplacian center gauge [16]. The proposal was to a large extent inspired by the 
Laplacian Landau [17], Laplacian abeUan [18], and Laplacian center [19] gauges. 
The idea is the following: 

To find the "best fit" to a lattice configuration by a thin center vortex configu- 
ration one looks for a matrix M{x) maximizing the expression: 

fKM = [M^{x)UAi,ix)M{x+fi)] , (13) 

x,ti 

with a constraint that M{x) should be an S0(3) matrix at any site x: 

M^(x) ■ M{x) = 1, detM(x) = 1 . (14) 
We soften the orthogonality constraint by demanding it only "on average": 

(M^ • M ) = — Y.M'^ W • M{x) = 1. (15) 

It is convenient to write the columns of M{x) as a set of 3-vectors: fa{x) = 
Mab{x). The optimal M{x), maximizing 1{m with the constraint, is determined by 
the three lowest eigenvectors /j" (x): 

'Dij{x,y)f'j{y)=Kfn^) (16) 
of the covariant adjoint Laplacian operator !D,y(x,>'): 

"DiMj) = ~ pAti{x)\ij^y,x+(. - [UA^,{x-fl)]jihy,x-(^ ■ (17) 

The resulting real matrix field M{x) has further to be mapped onto an SO(3)- 
valued field gA{x)- A naive map (which could also be called Laplacian adjoint 
Landau gauge) amounts to choosing gA{x) closest to M{x). Such a map is well 
known in matrix theory and is called polar decomposition. 
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OR ... Over-Relaxation 



Figure 1. Scheme of the Laplacian mapping of M{x) to a nearby S0(3) matrix field gA{^)- 

A better procedure, in our opinion, is the Laplacian map, that leads to direct 
Laplacian center gauge. We try to locate gAi^) as close to M{x) local maximum 
of the MCG (constrained) maximization problem. To achieve this, we first make 
the naive map (polar decomposition), then use the usual quenched maximization 
(overrelaxation) to relax to the nearest (or at least nearby) maximum of the MCG 
fixing condition. This procedure is illustrated schematically in Fig. 1 . 

To test the new procedure, we have recalculated the vortex observables intro- 
duced in our previous work (cf. Refs. [1, 2]), with P-vortices located via center 
projection after fixing the lattice to the new direct Laplacian center gauge. The 
results are summarized on the following page. 

The quantities of the most immediate interest are the center-projected Creutz 
ratios. Our data for the range of couplings P = 2.2 — 2.5 is displayed on a logarith- 
mic plot in Fig. 2a. In general Xcp{R,R) deviates from the full asymptotic string 
tension by less than 10%. 

As another way of displaying both center dominance and precocious linearity, 
we show, in Fig. 2b, the ratio 

Xphys {R,R)/Ophy, = x,p {R,R)/aLc,t ((3) (18) 

as a function of the distance in physical units Rphys = Ra{^) for all Xcp{R,R) data 
points taken in the range of couplings P = 2.3 — 2.5. Again we see that the center- 
projected Creutz ratios and asymptotic string tension are in good agreement (de- 
viation < 10%), and there is very little variation in the Creutz ratios with distance. 
We should probably stress in this context the significance of precocious linearity: 
it implies that center-projected degrees of freedom have isolated the long-range 
physics, and are not mixed up with ultraviolet fluctuations. 

Other encouraging results from MCG are recovered in the new gauge as well. 
Figure 2c shows the P-vortex density vs. (3 in a logarithmic plot. The density 
scales according to the asymptotic freedom formula with the slope corresponding 
to a quantity that behaves like a surface density. The slope for pointlike objects 
(hke instantons), or linehke objects (like monopoles) would be quite different. 
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Direct Laplacian Center Gauge 
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Figure 2. Results from direct Laplacian center gauge: 

(a) Combined data, at P = 2.2 — 2.5, for center-projected Creutz ratios obtained after direct Lapla- 
cian center gauge fixing. Horizontal bands indicate the asymptotic string tensions on the unprojected 
lattice, with the corresponding errorbars, taken from [20]. 

(b) The ratio of projected Creutz ratios to the full asymptotic string tension, as a function of loop 
extension in fermis. The data is taken from XcpiR,R) at a variety of coupUngs and lattice sizes. 

(c) Evidence of asymptotic scaling of the P-vortex surface density. The solid line is the asymptotic 

freedom prediction with \/ p /6A^ = 50. 

(d) Ratio of one- and two-vortex to zero-vortex Wilson loops Wi^2{C) /Wq{C) vs. loop area, at 
P = 2.3 on a 16'^ lattice. 

(e) Creutz ratios on the modified lattice, with vortices removed, at P = 2.3. 
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Figure 2d presents the data on vortex-limited Wilson loops. One can clearly 
see the expected trend, see Eq. (5), for large enough loops. Figure 2e shows that 
removal of center vortices causes the asymptotic string tension to vanish. 

6. How Does DLCG Differ from Laplacian Center Gauge? 

The first step of direct Laplacian center gauge fixing is similar to the Laplacian 
center gauge proposed by de Forcrand and collaborators [19]. Instead of using 
the three lowest eigenvectors of the covariant adjoint Laplacian operator and the 
naive map (or polar decomposition, see above), de Forcrand et al. build on the two 
lowest eigenvectors only. The gauge is fixed by g(x) that 

1. makes the lowest lying eigenvector to point in the third color direction (U(l) 
invariance still remains), and 

2. rotates the second lowest eigenvector into (say) the first color direction. 

There is an ambiguity in the procedure when the first and second vectors are 
coUinear, and such ambiguities should define positions of center vortices. 

LCG has its virtues and vices. It is unique (apart from eventual true Gribov 
copies) and shows center dominance after center projection. On the other hand, 
center dominance is seen only for very large distances, and there is not a good 
separation between confinement and short-range physics: there is no precocious 
linearity, there are too many vortices, vortex density does not scale. Moreover, 
identification of vortices via gauge fixing ambiguities fails for simplest configura- 
tions (like a pair of thin vortices put on the lattice by hand [12]), and is practically 
impossible in Monte-Carlo generated configurations. Center projection is neces- 
sary. 

To improve on these problems, Langfeld et al. [21] proposed to follow the 
LCG procedure of de Forcrand et al. by (over-)relaxation to MCG. This, in anal- 
ogy with DLCG, could be called indirect Laplacian center gauge. ^ 

The question is whether results from DLCG and ILCG differ considerably, and 
whether there is any correlation between vortex locations in those two gauges. 
Figure 3a shows the projected Creutz ratios at P = 2.4 in ILCG; for compari- 
son we also display the corresponding data from DLCG. It seems that the center 
dominance properties are somewhat better in DLCG than in ILCG, though the dif- 
ference is not great. The reason for this is quite easy to explain: both procedures 
seem to locate the same physical vortices. 

The simplest way to test the last statement is the following: For a given lattice 
{Ufji{x)} let {Z'i^{x)} be the lattice obtained by center projection in DLCG, while 
{Z'lJx)} be the corresponding lattice in ILCG. Denoting by Z'(C) and Z"(C) the 

^LCG involves first fixing to Laplacian abelian gauge, then further reducing the residual sym- 
metry from U(l) to Z2, in which it is reminiscent of indirect maximal center gauge of Ref. [1]. 



9 
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Figure 3. Comparison of the direct and indirect Laplacian center gauges: 

(a) Center projected Creutz ratios. (The asymptotic string tension is shown by the horizontal band.) 

(b) Creutz ratios XpwdiM^f^) calculated from "product Wilson loops", Eq. (19). 



Wilson loops in these two projected lattices, we construct the "product" loops 

Zprod{C) = {Z'{C)Z"{C)) (19) 

and from their expectation values the corresponding Creutz ratios Xprod{I,J)- The 
expectation is that if the two projected lattices were perfectly correlated 

Xprod{R,R) =0 (perfect correlation), (20) 

whereas in case of zero correlation 

Xprod{R:R)=y:!cp{R,R)+l!cp{R,R) (no correlation). (21) 

It is evident from Fig. 3b that center-projected loops in both gauges are not 
well correlated at short distances, but become correlated at large distances. The 
interpretation, we believe, is straightforward: The P-vortices in each projected 
lattice do not coincide, but in most cases are located within the same (thick) center 
vortices on an unprojected lattice. This accounts for the strong correlation on large 
distance scales, larger than the typical size of vortex cores. Similar correlations 
exist also with projected lattices in MCG (with gauge fixing via overrelaxation). 

7. Summary 

1. Center dominance exists in various gauges. The maximal center gauge has 
an appealing "best-fit" interpretation, but the successes of the approach have 
been overshadowed by the problem of Gribov copies. 

2. We have proposed a new gauge, direct Laplacian center gauge, that combines 
fixing to adjoint Laplacian Landau gauge with the usual overrelaxation. The 
first step of the procedure is unique, in the second step no strong gauge- 
copy dependence appears. This procedure can be interpreted as a "best fit" 
softened at vortex cores. 
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3. All features known from MCG are reproduced in direct LCG: center domi- 
nance, precocious linearity, scaling of the vortex density, etc. 

4. Similar results follow from center projection in Laplacian center gauge after 
overrelaxation (indirect LCG). The reason is that vortex locations in pro- 
jected lattices in direct and indirect LCG are quite strongly correlated. 
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